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We present our latest study of Lambda-Nucleon (AAO interaction by using lattice QCD, following 
up on our report at LATTICE 2008. We have calculated not only the scattering lengths but also the 
central and tensor potentials, which are obtained from the Bethe-Salpeter (BS) amplitude mea- 
sured in lattice QCD. For these calculations, we employ two different types of gauge configura- 
tions: (i) 2+1 flavor full QCD configurations generated by the PACS-CS collaboration at J3 = 1.9 
(a = 0.0907(13) fm) on a 32 3 x 64 lattice, whose spatial volume is (2.90 fm) 3 , with the quark 
masses corresponding to (m n ,m K ) « (301,592), (414,637), (570,724) and (699,787) (in units 
of MeV). (ii) Quenched QCD configurations at j3 = 5.7 ( a = 0.1416(9) fm) on a 32 3 x 48 lattice, 
whose spatial volume is (4.5 fm) 3 , with the quark masses corresponding to (m z ,mK) ~ (512,606), 
(464,586) and (407,565). The following qualitative features are found: The Ap potential has a 
relatively strong (weak) repulsive core in the l So ( 3 Si) channel at short distance, while the poten- 
tial has slight attractive region at medium distance. The tensor potential is found to be weaker 
than the AW case. These results hold in both full and quenched QCD. The energy of the ground 
state on the finite lattice volume is calculated: In both spin channels, energy shift due to the finite 
volume, from which we extract the scattering length via the Luscher's formula, is found to be 
negative at all values of quark masses, suggesting that the AN interaction is attractive. We have 
also discussed the quark mass dependences of the potentials and the scattering lengths. 
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1. Introduction 

The hyperon-nucleon (YN) and the hyperon-hyperon (YY) interactions are the keys to explore 
the strange nuclear systems, in which hyperons (or strange quarks) are embedded in normal nuclear 
systems as "impurities" [ID]. For example, spectroscopic studies of the A and £ hypernuclei are 
carried out experimentaly and theoretically, which leads to a qualitative conclusion that the A- 
nucleus interaction is attractive while the Z-nucleus interaction is repulsive. Such information is 
useful to study the composition of hyperonic matter inside the neutron stars [Q] : the A particle 
instead of Z~ would be the first strange baryon to appear in the core of the neutron stars. The 
EN interaction is also interesting and important, in order to explore the existence of E hypernuclei 
and the dense hyperonic matter in neutron stars. Despite their importance, YN and YY interactions 
have still large uncertainties because direct YN and YY scattering experiments are either difficult 
or impossible due to the short life-time of hyperons. 

Under these circumstances, the lattice QCD would be a valuable theoretical tool to make a 
first-principle calculation of baryon-baryon interactions. Previously, scattering parameters based 
on the Liischer's formula have been reported for the AW system [Q, ^J and for the YN system |0, [7]]. 
Recently, a new approach to the AW interaction from the lattice QCD has been proposed [||]. In this 
approach, the NN potential can be directly obtained from lattice QCD through the Bethe-Salpeter 
(BS) amplitude and the observables such as the phase shift and the binding energy can be calculated 
by using the resultant potential. See Refs. [^ [l(J 11, [l2|, |J, U for the recent developments in 
various aspects along the line of this approach. 



The purpose of this report is to present our recent calculation of the AW potentials together with 
the scattering parameters by using the full and quenched QCD gauge configurations. A preliminary 



result was already reported at LATTICE 2008 [15]. This is the latest version of the report which 



includes several new efforts: (i) Not only the central potential but also the tensor potential are 
calculated, (ii) quark mass dependences of the potentials and the scattering lengths are studied, and 
(iii) statistics of the lattice are increased from the previous report. 

2. Formulation 

The basic formulation has already been given in Refs. [|| |TT|, 16, 17]. (See also Refs. [|3, 18].) 



and a recent comprehensive accounts for the lattice AW potential is found in [ 14]. The formulation 



for the YN interaction is reported in [ 15 ], so that we mention only the basic equations briefly below. 



We start from an effective Schrodinger equation for the equal-time BS wave function: 

-J-V 2 <Kr)+ Ju{r,r f )${T , )d>T' = E<j)(r). (2.1) 

Here /x = mj<jntj / {m\+m^) and E = k 2 /(2ju) are the reduced mass of the AW system and the non- 
relativistic energy in the center-of-mass frame, respectively. We consider the low-energy scattering 
state so that the nonlocal potential can be rewritten by the derivative or velocity expansion Jl^], 
JJir^r 1 ) = VAN(r,V)8(r — f). The general expression of the potential Van is known to be [ |20| ] 

Van = Vo(r)+V a (r)(aA-ON) + VT(r)S l2 + V L s(r)(L-S + ) + V ALS (r)(L-S-) + 0(V 2 ). (2.2) 
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Here S12 = 3(cTa • n){c^ ■n) — <$a • &N i s the tensor operator with n = r/\r\, S± = (5)v ± Oa)/2 are 
symmetric (+) and antisymmetric (— ) spin operators, L = —if x V is the orbital angular momentum 
operator. Vq&j are the leading order (LO) potentials while Vls^ls are the next-to-leading-order 
(NLO) potentials in the velocity expansion. The procedure to derive the LO potentials is as follows. 
For the spin triplet state, the noncentral forces (e.g. the tensor force) mix different partial waves 
such as S- and D- waves. Namely, the spin-triplet wave function (j>(r;J = 1) comprises the S- and 
the D-wave components, which can be extracted from the lattice wave function <§> a p{r; J = 1) such 



that (See Q15[ ] for the method to obtain (j) a p(r,J) in lattice QCD.) 



(2.3) 



f a/5 (r; 3 Si) = ^<p a p(r,J = 1) = £E# eO #0 a /j(r,7 = 1), 
\ <j> a p(r; 3 Dj) = £<j> ap (r,J = 1) = (1 - &)* a p(nJ = 1). 

Therefore, the effective Schrodinger equation with the LO potentials becomes: 

I 2 \ x {-^ 2 + V °(r)+V a (r)(0 A -0 N )+V T (r)S l2 ^(r)=l ^\xE(f>(r) (2.4) 

The r dependence of the central and the tensor potentials, Vc(r;J = 0) = Vb(r) — 3V a (r) for 7 = 0, 
Vc (r; J = 1 ) = Vb (r) + V ff (r) , and Vr (r) for 7 = 1 , are determined once we obtain the wave function, 
the total energy and the reduced mass in lattice QCD. 

3. Numerical simulations 

3.1 2+1 flavor QCD 

Main results in this report are obtained by using the 2+ 1 flavor full QCD gauge configurations 
generated by PACS-CS collaboration [ pl| ] with the RG-improved Iwasaki gauge action and the non- 
perturbatively 0(a)-improved Wilson quark action at j3 = 1 .9 on a 32 3 x 64 lattice. The lattice spac- 
ing at the physical quark masses has been estimated as a = 0.0907(13) fm [p|]. So far, we have em- 
ployed four values of the hopping parameter for light quarks, k iic i = 0.13700,0.13727,0.13754,0.13770, 
while the one for the strange quark is fixed to K s = 0.13640. Several light hadron masses obtained 
in the present calculation are shown in Table ||. To calculate the BS wave function, the wall source 
is placed at the time-slice to with the Coulomb gauge fixing, and the Dirichlet boundary condition 
is imposed in the temporal direction at the time-slice t — to = 32. In order to improve the statistics, 
multiple sources at f = 8n with n = 0, 1,2, ■ • • ,8 are employed on each gauge configuration. The 
results are obtained with AUrf = 609,481,568,422 for K ud = 0.13700,0.13727,0.13754,0.13770, 
respectively, where Af con f is the number of the gauge configurations. 

3.2 Quenched QCD with larger spatial volume 

In quenched QCD calculation, we employ the plaquette gauge action and the Wilson quark 
action at j8 = 5.7 on a 32 3 x 48 lattice. The periodic boundary condition is imposed for quarks in 
the spatial direction. The wall source is placed at to = with the Coulomb gauge fixing and the 
Dirichlet boundary condition is imposed at t = 24 in the temporal direction. The lattice spacing 
at the physical point is determined as a = 0.1416(9) fm (1/a = 1.393(9) GeV) from m p = 770 
MeV. The spatial lattice volume is (4.5fm) 3 . The hopping parameter for the strange quark mass 
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is given by K s = 0.16432(6) from tuk = 494 MeV. Three values of the hopping parameter, k U( i = 
0.1665,0.1670,0.1675, are employed for the light quark mass. Table |l| also lists the light hadron 
masses calculated in quenched QCD. The results in quenched QCD are obtained with iV con f « 1000. 

4. Numerical results 



4.1 2+1 flavor QCD 

Figures [T] and ^| show the AN potentials obtained from 2+1 flavor QCD calculation as a func- 
tion of r. The central (Vc(J = 1)) and the tensor (Vj) potentials in the 3 S\ — 3 D\ channels are given 
in Fig. [T] while the central potential in the 'So channel (Vc(J = 0)) is given in Fig. [| We also show 
the central potential multiplied by volume factor (r V c {r)) in the left panel in addition to the nor- 
mal V (r) given in the right panel, in order to compare the strength of the repulsive force between 
two quark masses. These figures contain results with im n ^m K ) w (699,787) and (414,637) MeV, 
which are obtained at t — to = 13 and 10, respectively. These time-slices are chosen so that the 
ground state saturation is achieved. 

As can be seen in both figures, the attractive well of the central potential moves to outer region 
as the u, d quark mass decreases while the depth of these attractive pockets do not change so much. 



The present results show that the tensor force is weaker than the AW case []17[], and the quark mass 
dependence of the tensor force seems to be small. Both of the repulsive and attractive parts increase 
in magnitude as the u,d quark mass decreases. 

For m n ss 700 MeV, the central potentials reach Vc — > at the radial distance r ~ 1 .3 fm, which 
is smaller than the half of the physical lattice length (aL/2 « 1.45 fm). Therefore the Liischer's 
formula can be applied to extract the scattering phase shift, which will be discussed in the latter 
subsection. For m n « 400 MeV, on the other hand, the interaction range of the Vc, which is about 
1 .4 fm, almost reaches to the half of the lattice. Therefore we must be very careful to extract the 
scattering phase shift at this or lighter quark masses from the Liischer's formula, though no sign 
of the violation against the Liischer's condition was observed within errors for the effective central 
potential even at m % ss 300 MeV. (See Fig. 1 in the previous report [15].) Calculation on larger 
spatial volume will be needed to correctly extract the scattering phase shift at m % w 300 MeV. 



Kud 



m Tl 



m,, 



m K 



m K * 



m N 



m A 



mz 



m- 



2+1 flavor QCD by PACS-CS with K s = 0.13640 

0.13700 699.4(4) 1108(3) 786.8(4) 1159(2) 1572(6) 

0.13727 567.9(6) 1000(4) 723.7(7) 1081(3) 1396(6) 

0.13754 413.6(6) 902(3) 636.6(4) 1026(3) 1221(7) 

0.13770 301(3) 845(10) 592(1) 980(6) 1079(12) 



1632(4) 1650(5) 1701(4) 

1491(4) 1519(5) 1599(4) 

1349(4) 1406(8) 1505(4) 

1248(15) 1308(13) 1432(7) 



quenched QCD with j3 = 5.7, K s = 0.1643 
0.1665 511.8(5) 862(3) 605.8(5) 
0.1670 463.6(6) 842(4) 586.3(5) 
0.1675 407(1) 820(3) 564.9(5) 



Exp. 



898(1) 1297(6) 1344(6) 1375(5) 1416(3) 
895(2) 1250(9) 1314(9) 1351(6) 1404(4) 
886(3) 1205(13) 1269(9) 1326(9) 1383(5) 



135 



770 



494 



892 



940 



1116 



1190 



1320 



Table 1: Hadron masses in the unit of MeV. 
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Figure 1: The central and the tensor potentials in S\ — Dj channel in 2 + 1 flavor QCD as a function of r 
at m% ~ 414 MeV (red and blue) and 699 MeV (green and magenta). 
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Figure 2: The central potential in ^Sq channel in 2 + 1 flavor QCD as a function of r at m K ~ 414 MeV (red) 
and 699 MeV (green). 



4.2 Quenched QCD 

Figures | and |] show the AN potentials with (m K ,m K ) ~ (512,606) and (407,565) MeV, in 
quenched QCD at the time-slice t — to = 7. Fig. [3] shows the Vc(J = 1) and Vj, while the Fig. ^ 
shows the Vc(J = 0). We find that the qualitative behaviors of the AN potential in quenched QCD 
are more or less similar to those in full QCD in both 7 = 1 and channels: Namely, the attractive 
pocket of the central potential moves to longer distance region as the quark mass decreases, and 
the quark mass dependence of the tensor potential seems to be small. 



4.3 Scattering lengths 

Figure B| shows the scattering lengths as a function of m\, which are calculated through the 
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Figure 3: The central and the tensor potentials in 3 5i — 3 D\ channel in quenched QCD at m n ~ 407 MeV 
(red and blue) and 512 MeV (green and magenta). 



> 

CD 



20 


* 




v t * 

© v * 


10 



- o V 
• r 

<x> 

•I 





> 





50 
40 

30 - 
20 
10 
_ 

^ -10 
-20 
-30 
-40 
-50 



^=407 MeV 
m„=512MeV 



W« 



All 




0.0 0.5 

r(fm) 



0.5 



1.0 
r(fm) 



1.5 



2.0 



Figure 4: The central potential in 'So channel in quenched QCD at m K ~ 407 MeV (red) and 512 MeV 
(green). 



Liischer's formula [^ [T8|] 

A:cot6b(*) = ^Zoo(l; (^/(2?r)) 2 ) = l/a + O(£ 2 ), 
with Z 0Q (s;q 2 ) = ^Z HeZ3 (n 2 -q 2 )- s (Re ^> 3/2), 



(4.1) 



A- 



where Zoo(l;g 2 ) is obtained by the analytic continuation in s. The energy, E = £j, on the finite 
lattice volume is determined by fitting the asymptotic region of the BS wave function in terms of 
the Green's function 



1 



G(^) = 73l 



1 



jpr 



L3 ^P 2 ~ k2 



r~ \p,p = n — ,n£Z 3 }, 



(4.2) 



which is the solution to the Helmholtz equation (A + fc 2 )G(r,& 2 ) = — 5j,(r) with 5/,(r) being the 
periodic delta function^, p| ]. 
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As is seen in the Fig. |5], the scattering lengths are almost constant for larger u,d quark mass 
corresponding to m % > 560 MeV. On the other hand, for lighter u,d quark mass region that 400 
MeV < m K < 500 MeV, the present result seems to show that the scattering lengths increase as the 
u,d quark mass decreases. The present values are still much smaller than the empirical scattering 
lengths of AN, oq ~ 1.5 — 2.5 fm, estimated from the measurement of the AN total cross section 
and the theoretical studies of A-hypernuclei. The scattering lengths calculated in quenched QCD 
are qualitatively similar to those in 2 + 1 flavor QCD. As is discussed in the former subsection, we 
will need larger spatial volume in full QCD to extract the scattering lengths reliably at m K < 300 
MeV. 
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Figure 5: Scattering lengths of AN interaction as a function of m n . 



5. Summary 

We have calculated the central and tensor parts of AN potentials as well as the scattering 
lengths in lattice QCD at several values of the u,d quark mass corresponding to m n ~ 300 — 700 
MeV. The central potentials indicate that the interaction range becomes larger while the depth of the 
attractive well hardly changes as the u,d quark mass decreases. On the other hand the tensor force 
has relatively a weak quark mass dependence. The present result of the scattering lengths shows 
that the AN interaction is attractive and becomes stronger as the u,d quark mass decreases. The 
calculation in 2 + 1 full QCD for the AN system with larger volume and smaller lattice spacing at 
physical quark mass is highly desirable for definite conclusions on the nature of the AN interaction. 
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